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EVALUATION 


This  report  is  another  building  block  in  the  solution 
for  the  transient  response  of  either  a  metallic  or  dielectric 
targets.  These  solutions  will  provide  data  over  the  entire 
frequency  spectrum  and  will  reduce  the  need  for  extensive 
measurement  programs.  The  computation  of  the  response  of  a 
Tomahawk  like  target  will  be  useful  in  cruise  missile  detection 
and  identification  systems.  This  effort  has  application  to 
TPp^RlC  Surveillance  Sensor  Techniques. 
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SECTION  I 
INTRODUCTION 


This  is  the  final  report  under  Contract  No.  F30602-79-C-0060. 

The  main  result  is  the  application  of  the  space-time  integral  equation 
approach  to  the  dielectric  scattering  problem.  The  response  of  dielectrics 
is  of  recent  renewed  interest  due  to  the  potential  use  of  radar  absorbing 
material  (RAM)  on  radar  targets.  The  classical  frequency  domain  solution 
was  given  for  the  sphere  by  Mie  in  1906  (see,  for  example  ref.  Cl]) .  Modern 
computing  power  makes  short  work  of  such  a  solution;  see  for  example  the 
extension  of  the  Mie  result  to  layered  spherical  shells  in  ref.  [2].  In 
recent  years,  numerical  techniques  have  been  brought  out  to  free  the  problem 
from  the  constraints  of  the  separable  coordinate  system  imposed  by  the  clas¬ 
sical  solution.  Most  of  these  techniques  are  applied  in  the  frequency 
domain:  Mautz  and  Harrington  C 3 1  investigated  a  frequency  domain  solution 
using  surface  integral  equations  on  solids  with  axial  symmetry.  A  time 
domain  solution  is  presented  in  [43  for  complex  dielectrics,  with  simplified 
geometries.  The  time  domain  integral  equation  approach  has  been  pursued  at 
the  Sperry  Research  Center  for  a  number  of  years  [5,6].  Results  have  been 
obtained  for  conducting  solids,  conducting  plates  and  compound  targets. 

These  computations  have  been  verified  by  measurement  on  the  time  domain 
scattering  range.  In  the  present  effort  these  techniques  are  directly 
extended  to  dielectric  solids. 

In  Section  2  the  theory  of  the  space-time  integral  equation  approach 
is  presented.  The  numerical  implementation  is  detailed  and  results 
are  given  for  several  shapes.  The  response  of  the  sphere,  in  comparison 
with  the  classical  solution,  is  used  as  one  aspect  of  the  verification  of 
the  technique. 

In  Section  3  the  technique  and  solution  is  given  for  scattering  from 
a  complex  target  consisting  of  wings  on  a  fuselage.  A  model  closely 
resembling  a  Tomahawk  missile  is  used. 


The  calculations  are  verified  by  measurement  on  the  scattering  range. 
The  scattering  range  and  results  of  response  measurements  on  a  variety  of 
targets  are  described  in  Section  4. 


SECTION  II 

SPACE  TIME  INTEGRAL  EQUATION  SOLUTION  FOR  SCATTERING  FROM  DIELECTRIC  SOLIDS 


2 . 1  THEORY 


Consider  a  closed  surface  S  separating  a  uniform  exterior  (e,  y,  c)  from- 
a  uniform  dielectric  interior  with  =  e^e,  y^  =  y^y,  and  =  c/y^e^y^  as 
shown  in  Figure  2-1.  In  the  time  domain,  the  exterior  electromagnetic  field 
can  be  written  in  terms  of  its  retarded  values  on  S,  a  form  known  as  the 
Kirchhoff  integral  representation. 


H(r,t) 


*£/  j£fi'  x  <fi'  >■')  x6)j  ds' 


(2-1) 


E(r,t)  =  E  (r,t)  + 


-  ~  n'  x  +  i(( n’  •  E')R  +  (n'  xE’)  xR  )|  dS' 


(2-2) 


where 

E'  means  E(r’,t),  etc. 
T  =  t  -  R/c 


R  is  the  magnitude  of  R  =  r  -  r’ 


R  is  the  unit  direction  of  R 


L  = 


_1_  _9_ 
Rc  3t 


n'  is  the  outward  normal  at  r’ 


These  equations  are  derived  using  the  Green's  theorem  with  the  time  domain 
Green's  function  G(r,r'|t,x)  =  1/R  6(t  -  R/fc  -T)  and  making  use  only  of  the 
source-free  Maxwell's  equations 
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FIG.  2-1  Scattering  problem  geometry. 
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A  readable  derivation  in  the  frequency  domain  is  given  by  Jackson  [7],  The 
external  sources  are  represented  by  H  and  E  (the  incident  field) .  A  similar 
pair  of  integral  equations  can  be  written  for  the  interior  fields  and  E  . 
But  inside  there  is  no  incident  field  and  the  sign  of  the  integral  is  changed 
as  n  remains  pointing  outward: 


/  £  3e  11 

+1,( («' •»(•)*  +  («•*»!’)  »S) 


l  /  ■p  *  * 

---sr  /  -is'  » e) 


dS1 

(2-3) 

dS ' 

(2-4) 


for  r  inside. 


where 


E|'  means  E 


\) 


Ti  "  *  -  R/Ci 


-  J_  _L_  jl 

1  “  R2  +  RC1  3t 


The  double  prime  is  used  to  emphasize  that  the  retardation  is  at  speed  c^. 
These  equations  permit  the  calculation  of  the  field,  in  particular  the  far- 
scattered  field,  from  the  surface  quantities  once  these  are  known.  Our 
problem  is  to  solve  the  equations  on  the  surface. 


! 

I 

1 

i 


! 
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For  r  on  S  the  term  with  the  operator  L  in  each  of  the  integrals  is 
singular.  This  singularity  can  be  removed  in  the  following  way:  consider 
equation  (2-1) .  We  note  that  the  integral  takes  on  the  same  value,  if  the 

— V 

singular  point  is  excluded,  whether  r  is  just  inside  or  just  outside  of  S. 
However,  the  singular  part  changes  sign  as  r  passes  through  S.  We  now  re¬ 
quire  that  equation  (2-1)  have  value  zero  evaluated  inside  S.  That  is, 


■+  -*  +  -H  ->•  -*■ 

H  (r  e  S  )  =  H  +  I  +1 
NS  S 


0  =  H  (r  e  S  )  =  H  +  I  -I 

NS  S 


Adding,  we  get 


H(r  e  S)  =  2H  +21 


Here,  I  is  the  integral  with  the  nonsingular  point  removed;  usually  written 
NS 

as  /  . . .  dS’.  A  similar  argument  is  made  with  equation  (2-3) ,  where  now 
"►  “►  + 

H^r  e  S  )  =  0.  This  is  no  more  than  a  restatement  of  the  properties  of  a 
Green's  function  solution  to  a  boundary  value  problem.  The  discontinuity 
at  the  boundary  S  is  then  required  to  be  expressed  by  the  known  physical 
boundary  conditions.  These  boundary  conditions  are 


nxH^  =  nxH  =  J 


E^xn  =  Exfi=M 


y(n,H1)=n*H=H 
r  \  1/  n 

Er(fi  •  Ex)  =  fi  •  E  E  En 


The  first  of  these  boundary  conditions  also  serve  to  define  the  tangential 

electric  and  magnetic  surface  currents.  (Also  one  could  define  the  surface 
“►  A  A 

charqes  o  =  E  •  n  and  a  =  H  •  n,  but  here  we  will  prefer  to  write  E  and 
3  m  n 

H  .)  For  convenience,  we  will  also  normalize  the  exterior  constants 
n 

e  =  y  si  c  =  1  so  that 
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and 


(2-10) 


/Si.  E  [Jl  k),  e 

/  cl  r\/  y  c  /  r 


ei  ci  rv  e  c/  r 


With  these  simplifications  and  application  of  equation  (2-5)  and  the  boundary 
conditions,  we  obtain  the  following  set  of  equations  for  r  on  S: 


J(r,t)  =  n  x  H 


•  2V  -J  1  /  (  Er  9M ’’  /Hn 

r,t)  =  n  x  H2  =  -  n  x  —  j-  j-  —  —  +  ^ 


R  +  J"  x  R  | !  dS '  (2-12) 


M(r,t)  =  -  n  x  E  =  2 


->  -+ 

M(r,t)  = 


nxE  =  2M1-nx  f  |‘  £  +  L(En  *  "  x  *)  J  ds'  (2-13) 

1  11  \  "  "  11  H  /  j'  TT  ¥  +  ‘1(77  *  -  ""  11  ■’>)}  ds'  (2-14> 


The  above  are  for  the  tangential  components.  For  the  normal  components 


we  can  write: 


H  =  2H1  +  n  •  f  j 
n  n  2TT  J  I 


2*  /  {*•• 


. *  as  in  2-11  . . .  \  dS ' 


as  in  2-12  . . •  \  ds 1 


E  *  2E1  +  n  •  -i-  -f  j. 
n  n  2TT  J  \ 


. .  as  in  2-13  . . . (  dS ' 


(2-15) 


(2-16) 


(2-17) 


. .  as  in  2-14  . . . (  ds ' 


(2-18) 
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->i  _  A  -4x  ->  -►i  i  x 

where  J  =  n  x  H  (r,t)  and  similarly  for  M  ,  E  ,  M 

n  n 

Several  observations  are  in  order.  First,  (2-11)  through  (2-18) 

-4  “► 

are  too  many  equations.  The  field  vectors  E  and  H  are  not  independent  but 

related  through  the  Maxwell  equations.  Hence  equations  (2-11)  and  (2-12) 

for  the  surface  components  of  H  and  equation  (2-13)  for  the  normal  component 

of  H  are  sufficient  to  solve  the  problem.  Such  a  formulation  is  detailed  in 

Appendix  7.1.  This  approach  led  to  instabilities  in  the  numerical  solution 

and  was  abandoned.  Second,  it  is  known  that  only  the  tangential  components 
“►  “►  ->•  ->■ 

of  E  and  H,  that  is  J  and  M,  are  required  to  completely  specify  the  field. 
The  normal  components  can  be  eliminated  by  means  of  the  continuity  relations 

7  •  J  +  4:  (e  -  E  \  =  0  , 

s  dt  V  n  In/ 

(2-19) 

■  r\ 

V  ‘M  +  -5-/H  -h1)  =  o. 

s  dt  \  n  In/ 

The  operator  V  •  is  the  surface  divergence.  The  numerical  evaluation  of 

V  •  J  on  an  arbitrary  three  dimensional  surface  can  be  very  inaccurate, 
s 

For  that  reason,  the  explicit  determination  of  the  normal  components  (equa¬ 
tions  (2-15)  through  (2-18))  was  retained.  There  was,  at  the  outset,  some 
worry  that  the  direct  evaluation  of  and  E^,  without  application  of  the 
constraints  (2-19) ,  would  lead  to  instabilities  in  the  time  domain  solution. 
However,  the  present  formulation  appears  to  be  successful. 

Since  we  wish  to  retain  the  complete  formulation  of  equations 
(2-11)  through  (2-18) ,  we  need  to  combine  these  somehow  into  fewer  equations 
such  as  to  permit  the  unique  solution  for  the  six  scalar  components  of  J, 

M,  H^,  E^.  We  achieve  this,  using  a  general  technique,  similar  to  that 
described  by  Mautz  and  Harrington  C 3 1 .  Namely,  multiply  the  second  of  the 

-►  -f  4 

H  (or  J)  equations  by  g  and  add  to  the  first;  multiply  the  second  of  the  E 

4 

(or  M)  equations  by  a  and  add  to  the  first.  For  example,  from  (2-15) , 

(2-16)  we  obtain 
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Mautz  and  Harrington  showed  that  the  solution  obtained  is  unique  if  g,  J§ 

* 

satisfy  a  g  =  positive  real.  Their  formulation  was  in  the  frequency  domain. 
In  the  time  domain,  complex  a  means  a  =  3/3t.  The  above  requirement 

thus  becomes  that,  for 

2  •  °i +  “2  ii 
6  -  +  £  • 

the  conditions  8;Ja2  -  ct.^  =  0  and  >  a^2  guarantee  uniqueness.  It  is 

likely  that  other  non-zero  values  of  a,  jg  will  also  lead  to  unique  solutions. 
Here,  we  will  use  the  simpler  form  of  the  constants  a  -  and  jg  =  8^. 

The  solution  is  now  formally  complete.  To  solve  the  scattering  problem, 
we  first  solve  the  problem  for  all  time  on  the  surface.  This  last  can  be 
done  by  a  simple  marching  in  time  procedure  since  the  quantities  appearing 
inside  the  integrals  need  to  be  evaluated  at  the  retarded  times  T  and  x  . 

If  it  is  arranged  that  X  and  X  are  always  at  least  one  numerical  time  step 
earlier  than  the  current  time,  then  no  matrix  inversion  needs  to  be  performed. 
To  get  started,  the  incident  pulse  is  taken  as  the  time  and  bandlimited 
function 


The  process  is  started  at  a  time  sufficiently  early  for  the  above  to  be 
considered  zero. 

Having  completed  the  surface  calculation,  the  far  field  cam  then  be 
computed  quickly  for  the  backscatter  case  as  well  as  at  any  number  of 
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bistatic  angles,  by  a  limiting  form  of  equation  (2-1)  or  (2-2).  As  r  00 , 

the  scattered  field  is  given  by 


■*s 
r  H 


It  t  \  -  ±  i 

*  (  3M1 

3J‘  .  } 

4H  J 

|  3t 

+  IT  x  r  j 

dS' 


(2-21) 


where 


-+  ■>  -*• 

J*  -  J(r' ,  t') 

t*  *  t  ♦  r'  •  r  . 

The  far  field  time  t^  is  referenced  to  t  =  0  and  the  origin.  (That  is,  an 
impulse  reflected  from  the  origin  at  t  =  0  arrives  in  the  far  field  at 
tf  =  0.) 

-»s  A  3m 

Note  that  since  H  must  be  perpendicular  to  r,  the  integral  of  •  r 

over  the  closed  surface  S  must  vanish.  Hence  3M/3t  may  be  replaced  by 

(?  x  3M/3t)  x  r  in  the  equation.  This  observation  is  useful  when  determining 

the  relative  importance  of  different  regions  of  S  and  for  debugging  purposes. 

2 . 2  NUMERICAL  IMPLEMENTATION 

The  surface  s  is  divided  into  patches  of  approximately  equal  size,  with 

-► 

sample  points  at  their  centers.  Time  also  is  quantized.  The  values  of  J 
-► 

and  3J/3t,  etc.,  are  assumed  to  be  constant  over  a  patch.  In  evaluating  the 

/v 

integrals,  the  quantities  1/R  and  the  operators  L  R  are  also  assumed  to  be 

constant  oyer  the  patch,  except  for  the  "self-patch".  The  self-patch  (the 

■>  -► 

patch  containing  both  observation  and  integration  points,  r  =  r')  is  inte¬ 
grated  analytically.  The  result  is  removed  from  the  integral,  leaving  the 
"non-self"  integral. 

The  self- integrals  are  of  three  types.  To  evaluate  these  we  consider 
each  (curvilinear  rectangular)  self-patch  with  area  AS  as  a  circular  patch 
radius  Y-  Then  the  first  of  the  self  integrals  is  straightforward 
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_1_ 

2TT 


(2-22) 


The  others  require  the  expansion  of  ri  and  R  in  surface  coordinates.  This  is 
done  in  Appendix  7.3.  The  results  are 


with 


A 


n 


R  dS' 


e  H 
s  n 


(2-23) 


e 

s 


and 


n  x 


-L  / 
2tt  r 

•'as 


j' 


x  R  ds  =  e 


The  dyadic  e  is  no  more  than  a  convenient  short-hand  for  the  result,  which 

4  ► 

is  polarization  dependent.  The  magnitude  of  e  is 


I  hr  - 


(’V'v) 


(2-24) 


where  K  and  K  are  the  principal  curvatures,  and 
u  v 


e  *  e 


fa  3  -  a  a  I 

[u  u  V  V  I 


so  that 


e  •  J  =  e  a  J  -  e  a  J 
u  u  v  v 

Note  that  e  ■  0  for  a  sphere. 

The  standard  geometry  identifies  (a^,  av,  n)  as  (§,  ?)  for  the 

sphere.  For  a  body  of  revolution  about  the  z-axis,  au  becomes  the  tangent 
t  almost  parallel  to  9:  (t,  $,  n) .  Hence,  for  a  cylinder  of  radius  a, 
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we  have  on  the  cylinder  body,  =  0,  Ky  =  -  1/a,  so  that  e  =  J  (1/a)  >  0, 


and  e  =  -  (1/a)  <  0. 

s  4 


We  will  define  and  write  out  the  numerical  form  of  the  non-self  inte- 
grals  explicitly.  At  time  t^,  for  observation  point  r^: 
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(2-26) 


(2-27) 


(2-28) 


where 


M ' ,  etc . ,  are  evaluated  at  T  =  t .  -  R . . 
k  j  ik 


M."  ,  etc.,  are  evaluated  at  T,  =  t.  -  R.,/c  . 
k  1  3  ik  r 


At  each  time  step  t.  the  above  integrals  are  evaluated  for  each  observation 

point  i.  To  do  this,  it  is  required  to  keep  past  values  of  the  six  surface 

quantities  in  storage  for  all  times  back  to  R  /c  +  3At,  where  R  is  the 

max  r  max 

greatest  target  diameter,  and  At  is  the  time  step.  The  evaluations  at  T 
(or  T^)  are  achieved  by  polynomial  interpolation  over  values  at  the  five 
time  points  nearest  T.  (An  alternate  interpolation  method,  which  performs 
smoothing  on  the  five  points,  has  also  been  used  and  has  the  advantage  of 

suppressing  instabilities,  though  at  the  expense  of  a  somewhat  reduced 

» 

accuracy  in  the  final  result.) 

For  convenience  in  writing,  we  drop  the  subscripts  i,  j  in  the  following. 
The  numerical  forms  of  equations  (2-11)  and  (2-12)  are  thus 


-1 


(2-29) 


(««— ►  -4— ►  \  -►  _  3m  .  -► 

I  -  e  )  •  J  +  Ynxgr=2J  +  n  x  I 

(  I  +  e  )  •  J  -yernx^-=-nxl1 

As  described  in  Subsection  2.1,  we  will  multiply  the  second  of  these  by  Jg 
and  add,  obtaining 

(1  +  jg)  J  +y(l  -  £rg)  n  x  |^  =  2  J1  +  n  x  (i  -  Jg  ij  (2-30) 

Here  we  have  written  J  for  convenience  to  mean 

->■  ->• 

for  a  sphere  J  =  J;  in  any  case  this  is  a  small  correction  factor. 

The  pair  of  E  equations  (2-13)  and  (2-14)  will  result  in  a  similar 
form,  but  for  M  and  n  x  3j/3t.  To  solve  simultaneously,  we  thus  need  to 
write  a  numerical  form  for  the  time  derivative.  Here  we  choose  the  first 
order  form 


I  -  -  3(Vi))/4t  • 

Since  the  constants  y,  SS»  J§  will  be  such  that  the  3J/3t  and  3M/3t  terms  are 
relatively  small,  the  above  approximation  is  good  enough.  For  greater  accu¬ 
racy,  however,  one  could  use  the  second  order  form 

ff  -  (l.5  J(tj)  -  2  3(Vl)  *  .5  .  (2-33, 

We  now  multiply  (2-14)  by  a  and  add  to  (2-13) .  Also,  we  will  cross  multiply 
the  result  by  n,  in  order  to  obtain  the  same  unknowns.  (Note  that 
(n  x  M)  x  n  *  M,  since  M  is  a  surface  vector.)  Also  writing  the  time  deri¬ 
vative  numerically,  we  obtain  the  following  pair  of  equations  in  unknowns 
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J  and  n  x  M: 


2t 


I1  »  6)1  ♦  £  (i  -S*t)  fi  k5-  2  31  *  a  x(T  -gij  ‘ifi-tgS.,' 

( 1  ■  2  u^. )  J  +  (l+a)  nxJB-nx^M*-  n  x  |k  -  g  Kj  +  ^6  |l  -  <jUr) 


(2-3- 


The  right-hand  sides  of  these  are  just  combinations  of  the  integrals, 
already  expressed  in  (2-25)  through  (2-28) ,  with  the  incident  field  and  the 
previous  computed  values.  We  can  write  these  equations  as 


■»  m 

-*■ 

J 

V, 

1 

Q 

= 

. 

-¥ 

fi  X  M 

m  _ 

L  2  J 

and  solve  for  each  of  the  surface  components  u  and  v.  [in  the  present  com¬ 
puter  program,  the  vectors  are  stored  in  Cartesian  form,  so  that  the  equa¬ 
tions  are  solved  for  the  three  Cartesian  components.  This  introduces  a 
slight  complication,  since  J  as  well  as  appears  in  the  equation;  this  ef- 
feet  is  neglected  since  it  is  small.]  After  solving  for  J  and  M,  the  self¬ 
term  correction  dyadic  is  applied: 
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(2-36] 


with  a  similar  set  of  expressions  for  M  and  M  . 

u  v 


The  solution  for  the  normal  components  is  more  direct  as  (2-15)  and 


(2-18)  can  be  written  as 


where  e  is  defined  in  equation  (2-23) 
s 
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(2-37) 


Note  some  interesting  things  about  the  solution  (2-34) .  Suppose  we 
have  |$  =  1/e^  and  ot  =  1  /Vr#  then  the  cross  terms  disappear.  This  is  con¬ 
venient,  because  the  terms  which  have  y  in  them  are  the  result  of  the  self- 
patch  calculation.  They  are  large  and  inherently  inaccurate.  (The  self¬ 
patch  term  involving  e  is  small.)  The  equations  thus  reduce  to 


(2-38) 


That  is,  by  choice  of  the  weighting  factors  a  and  g,  we  were  able  to  elimi¬ 
nate  the  annoying  large  self-patch  contribution.  We  might  expect  this 
choice  of  a  and  £  to  give  the  best  results.  In  practice  it  was  found  that 
results  for  the  sphere  were  most  accurate  and  stable  for  a  =  1/y  =  1,  but 

]3  *  .7  1  /Er* 
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We  also  note  that  the  limit  e  -*■  00 ,  y  0  such  that  e  y  remains 

r  r  r  r 

finite  represents  the  conductor.  In  fact,  equations  (2-34)  do  reduce  to 
the  single  equation  used  in  [5,6]  and  Section  3  for  the  conducting  case: 

J=2JX+hxI.  (2-39) 

2.3  RESULTS  OF  THE  STIE  CALCULATION 

The  integral  equation  solution  was  tested  for  the  shapes  illustrated  in 

Figure  2-2,  namely  the  sphere  and  sphere-capped  cylinder,  both  with  =  3. 

For  the  sphere,  computational  results  could  be  compared  directly  with  the 

results  from  the  classical  solution.  This  comparison  is  given  in  Figure 

2-3.  The  incident  pulse  was  the  smoothed  impulse  of  equation  (2-20)  with 

a  =0.5.  As  can  be  seen,  the  agreement  is  nearly  exact.  It  is  interesting 
n 

to  compare  the  response  with  that  of  a  conducting  sphere,  shown  in  Figure 
(2-2b)  with  a  vertical  scale  reduced  by  a  factor  of  2. 

The  computed  response  of  the  sphere-capped  cylinder  is  shown  in  Figures 
2-4  and  2-5  for  TE  and  TM  polarizations  and  at  various  angles  of  incidence. 
The  second  bump  in  these  responses  is  interpreted  as  the  result  of  an 
internally  propagated  wave,  which  is  partially  reflected  at  the  backside  and 
reradiated  upon  re-emergence  from  the  front. 

In  Figures  2-6  and  2-7  are  shown  the  results  of  comparing  these  cal¬ 
culated  responses  with  those  measured  on  the  scattering  range  (section  4) . 
The  measured  responses  were  here  transformed  by  a  convolution  procedure  to 
obtain  the  response  to  the  standard  incident  pulse,  thus  permitting  direct 
comparison. 

The  parameters  of  the  sphere  were:  8  bands  of  approximately  equal 

sized  patches,  with  a  total  of  40  patches  on  the  hemisphere  (plane  symmetry 

was  utilized  requiring  the  computation  of  the  hemisphere  only) ;  radius  =  1; 

e  =  3.0,  y  =  1.0.  The  time  step  of  the  solution  was  At  =  .333,  while  the 
r  r 

solution  parameters  a  were  .22  and  1.0.  The  actual  values  of  6  and  a 
had  only  a  weak  effect  on  the  solution. 


J 
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(a)  SPHERE 


(b)  SPHERE-CAPPED  CYLINDER 

79-940 

FIG.  2-2  Test  targets  for  SHE  solution. 


FIG.  2-3  Smoothed  impulse  response  -  sphere: 
an  =  1/2,  er  =  3  (STIE). 


FIG.  2-7  Comparison  of  measured  and  calculated  responses  -  dielectric  sphere 
capped  cylinder  (er  =  3). 
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The  sphere  capped  cylinder  parameters  were  the  same,  except  for  the 
insertion  of  8  bands  of  8  patches  each  representing  the  cylinder  body. 


2.4  CLASSICAL  SOLUTION  FOR  DIELECTRIC  SPHERES 

The  eigenfunction  solution  for  the  sphere  is  given  by  Stratton  [1] 
among  others.  The  method  consists  of  expanding  the  incident  plane  wave, 
the  scattered  wave,  and  the  interior  wave  in  terms  of  spherical  wave  func¬ 
tions.  The  latter  form  a  complete  set  of  independent  solutions  to  the  wave 
equation.  The  boundary  conditions  are  then  applied  and  the  coefficients 
of  like  functions  equated.  We  collect  here  the  pertinent  equations.  The 
geometry  is  still  that  of  Figure  (2-1) ,  except  that  the  target  is  a  sphere. 

The  incident  plane  wave  from  the  positive  z  direction  can  be  written 
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The  spherical  wave  functions  are 
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In  the  preceding  rather  complicated  expression,  the  subscript  o  means  odd, 
e  means  even.  On  the  right-hand  side  the  upper  of  sin  4>  or  cos  4>  refers 
to  the  odd  form,  the  lower  to  the  even  form;  similarly  the  lower  of  the  ± 
sign  refers  to  the  even  form.  The  quantity  [...]*  means  3/3 (kr)  [...]. 

The  wave  number  k  •  =  w/c  =  w/ep"  . 

The  superscript  (1)  refers  to  the  appearance  of  the  first  kind  spherical 

Bessel  function,  j^Ckr).  The  first  kind  Bessel  function  is  finite  at  the 

origin  and  is  therefore  used  in  the  incident  wave  as  well  as  the  interior 

wave.  The  third  kind  spherical  Bessel  function  h^(kr)  appears  in  the 

n 

scattered  wave  expression,  since  it  is  outward  traveling  and  finite  at 

■+  ( 3 ) 

infinity;  its  appearance  in  the  wave  function  will  be  indicated  as  mQ^n ,  etc. 

The  function  P^  (cos  0)  is  the  Associated  Legendre  polynomial.  It  is  given 
n 


d  P  (cos  0) 

p  (cos  0)  =  sin  0  — — - — ~ 

n  d  (cos  0) 


where  the  Legendre  polynomial  is  found  recursively  by 


n  P  (x)  =  x  (2n  -  1)  P  (x)  -  (n  -  1)  P  (x)  (2-45) 

n  n-1  n-2 


the  scattered  field  is  written 
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In  the  interior,  we  have  similar  expressions  ‘for  E  and  H  (transmitted 
wave),  except  that  e,  p  ,  c,  k  are  replaced  by  e  ,  p  ,  c^,  k^;  the  coeffi¬ 
cients  are  a*"  and  bfc  and  the  eigenfunctions  contain  spherical  Bessel  func- 
n  n 

tions  of  the  first  kind.  We  now  evaluate  these  equations  at  the  boundary 
r  =  a  and  apply  the  boundary  conditions ,  namely 


r  X  \H  +  H  |  =  r  X  H 


(The  boundary  conditions  on  the  normal  components  are  noc  required.)  The 
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result  is  a  set  of  4  equations  in  the  four  unknowns  a  ,  b'5,  a  ,  b  ,  for 
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where  x  =  ka,  x^  =  k^a, 
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s  s 

The  above  are  solved  directly  for  a  and  b  , 

n  n 


the  coefficients  of 

the  scattered  field.  For  the  far  scattered  field  we  use  the  asymptotic  form¬ 
as  r  ■*■<*>.  For  instance,  at  4>  =  0 


•<3) 


rr  ’  =  -h(1)  (kr) 
eln  n 


dp 


.ikr 


d8 


-*■  1 


kr 


(-i) 


dP1 
n  n 


de 


.  ->(3) 
in,  =  l 
oln 


|kr  h^X) (kr)  J ' 


.ikr 


kr 


sin  9 


-*■  l 


kr 


(-i) 


n 


sin  0 


-24- 


For  the  other  polarization,  we  use  <J>  =  90  ;  then 
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Hence,  the  (normalized)  result  in  the  <f>  =  0  plane  (E-plane) 
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so  that  the  term  in  brackets  becomes  [a  -  b  ] 

n  n 


2.5  SOLUTION  FOR  TARGETS  WITH  DISPERSIVE  CONSTITUTIVE  PARAMETERS 

In  the  above  consideration,  the  constitutive  parameters  of  the  scat¬ 
tering  body  (y^,  G^)  were  taken  to  be  independent  of  frequency.  That  is, 
the  impulse  response  of  the  material  is  an  impulse.  For  many  dielectric 
materials  that  might  be  used  for  target  construction  this  assumption  is 
valid  over  the  range  of  operating  frequencies  that  may  be  encountered  in 
practice.  However,  there  also  exist  materials  and  situations  where  this 
assumption  of  frequency  independence  is  not  valj-d. 


2.5.1  Radar  Absorbing  Materials  on  Conducting  Surfaces 
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Radar  Absorbing  Materials  (RAM) ,  for  example,  typically  possess  con¬ 
stitutive  parameters  which  are  strongly  dependent  upon  the  frequency  of  the 


illuminating  waveform.  In  this  case  the  RAM  is  attached  to  the  skin  of  a 
conducting  target  and  its  effect  may  be  accounted  for  by  use  of  a  surface 
impedance  which  couples  the  tangential  E  and  H  field  on  the  surface  by: 


n  x  E  =  Z 


x  H)  x 


(2-53) 


where 


Zc  =  the  impulse  response  of  the  RAM  surface  impedance 
*  =  the  convolution  operation 

Note  that  in  this  case  there  again  is  no  penetration  of  the  fields  inside 

the  target  by  virtue  of  the  fact  that  it  has  a  conducting  skin.  These  two 

features  reduce  the  problem  to  the  solution  of  one  integral  equation  that 

is  derivable  from  (2-1)  by  use  of  the  impedance  boundary  condition  (2-53) 

A  — ► 

and  the  boundary  condition  that  n  •  H  =  0.  This  space-time  integral  equation 
is  given  by: 


n  x  H(r,t)  =  2n  x  H 
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T  =  t  -  R 


where 


Z  *  H  =  1  H(t)  Z  (t  - 1)  dx 

c  7  _  c 


The  solution  of  the  space-time  integral  equation  in  (2-54)  is  carried  out 
in  the  same  manner  as  has  been  used  for  the  case  of  perfectly  conducting 
targets.  The  problem  is  solved  numerically  for  the  tangential  surface 
field  n  x  H  by  marching  the  solution  on  in  time.  In  this  case,  however, 
a  second  term  appears  within  the  integral  on  the  right-hand  side  of  (2-54) . 
It  is  this  term  that  will  tend  to  reduce  the  tangential  field  setup,  and 
thus,  the  scattered  return.  The  value  of  this  term  results  from  the  con¬ 
volution  of  the  H  field  with  the  impulse  response  of  the  PAM.  This 
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convolution  can  be  computed  numerically  in  the  time  domain  since  it  uses 
only  previously  computed  values  of  H. 


2.5.2  Lossy  Materials 


Lossy  materials  may  be  modeled  by  allowing  an  electric  current  to  flow 
within  the  material.  In  this  model  Maxwell's  equations  are  written: 


r-t 

V  X  H  =  J 


„  3B 

v  x  E  =  -  aT 


(2-55) 


For  the  case  where  linearity  exists,  then  the  constitutive  relations  are 


D  =  £  *  E 


B  =  u  *  H 

-*  -* 
J  =  0  *  E 


(2-56) 


where  e,  u,  and  0  represent  the  impulse  response  of  the  constitutive 
parameters.  If  these  impulse  responses  are  impulses,  then  the  convolution 
becomes  a  simple  product.  Moreover,  the  problem  is  still  dispersive. 
Equation  (2-55)  becomes 


_  "►  9E  *+• 

VxH  =  E3t  +  0E 


3H 


V  x  E  =  -  v  t— 


3t 


(2-57) 


The  Fourier  transform  of  these  equations  give  their  frequency  domain  counter 
part  as 


V  x  H  =  (ju>  e  +  0)  E 

V  x  ?  =  -  ju  ij  S 


(2-58) 


/ 
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where 


H(r,t)  +-*■  HU,u>) 

E(r,t)  ■*-*•  E( r,w) 

and  an  equivalent  permittivity  e  may  be  written 

£  =  cjl  +  — — J  (2-59) 

V  3“ £ / 

The  time  domain  solution  of  this  problem  may  be  carried  out  in  several  ways. 

One  approach  would  be  to  solve  the  frequency  domain  integral  equation 
for  the  complex  e  at  the  frequencies  contained  within  the  spectrum  of  the 
incident  waveform  and  do  the  inverse  transform  of  the  weighted  result.  This 
approach  is  both  computer  time  and  memory  intensive. 

A  second  approach  would  be  to  solve  Maxwell's  equations  (2-57)  directly 
in  the  time  domain  throughout  the  volume  of  space  that  encloses  the  target 
and  its  neighborhood  using  an  extension  of  the  technique  first  proposed  by 
Yee  [9]  and  subsequently  by  Taflove  [10],  This  approach  is  computer  memory 
intensive  and  yields  only  approximate  results,  since  workers  to  date  have 
been  unable  to  account  exactly  for  the  effect  of  the  solution  volume  boundary 
on  the  numerical  results. 

A  Third  approach  would  be  to  develop  and  solve  the  space-time  integral 
equation.  The  frequency  domain  integral  equations  for  the  fields  inside 
the  body  are 
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E1  (r,u>) 


+  [^  +  i  (jt0/c)]  [("’  *®l) +  ("’ x  *l) x  R]j 


e”j“  R/c  ds- 


In  this  case 


1  =  .ATs  +  -£— \ 

c  /II  \  3  “  £1  / 


Consider  the  case  where  losses  are  low  so  that 


u>  e. 


<<  1 


and  then  (2-62)  becomes 


1  _  _1 

c  c 


1  (X  +  iu2El  ) 


where 


^-/vT 


Substitution  of  (2-63)  in  (2-60)  and  (2-61)  yields 


H1( r,w)  = 


’  *F  jf  j[i  (jwei +  °) x  ?i] 


•i.  +  i£ 

R2  R  2  /  C1 


1  +  i  i“ 

R  c. 


("'  ‘  *J  +  ("’  x  x  *] 


-ju>  R/c  “R\2  /  e 

e  1  e  '  '  ds 1 
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(2-61) 


(2-62) 


(2-63) 


(2-64) 


i 

J 


t  K 


The  space-time  integral  equations  are  now  obtained  by  taking  the  Fourier 
transform  of  (2-64)  and  (2-65) .  The  result  is 


where 


These  space-time  integral  equations  correspond  to  those  given  in  (2-3)  and 
(2-4)  for  the  lossless  target  case.  The  solution  is  carried  out  numerically 
in  the  same  manner  as  was  described  earlier  by  stepping  on  in  time.  It  is 
interesting  to  note  that  the  effect  of  losses  is  to  add  an  attenuation  factor 
e-aR  to  the  integrand  of  (2-66)  and  (2-67) .  There  also  appears  an  attenua¬ 
tion  term  in  the  L operator  which  is  proportional  to  the  conductivity  o. • 
This  solution  will  be  valid  for  the  case  where 


e 


»  O  E, 


(2-68) 


when  this  is  not  valid  then  the  direct  solution  of  Maxwell's  equations  in 
the  time  domain  appears  to  be  the  best  approach. 


2.5.3  Dispersive  Dielectric  Materials 

In  this  section  the  case  of  dielectric  materials  whose  permittivity 
varies  with  frequency  is  considered.  That  is 


e  =  e(w) 


so  that  in  the  time  domain 


D  =  e  *  E 


For  simplicity  the  magnetic  properties  are  constant  and  the  conductivity 
is  zero: 


y  =  constant 
a  =  0 


For  this  case  Maxwell's  equations  become 


V  x  H  -  e 


3E 

3t 


V  x  E 


(2-69) 
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where 


e  *  E  = 


E(x)  e (t  -  t)  dx 


(2-70) 


One  approach  is  to  solve  (2-69)  numerically  using  (2-70)  in  a  volume  of 
space  that  includes  the  target  and  its  neighborhood  by  marching  on  in  time. 
A  second  approach  is  to  solve  the  frequency  domain  integral  equations  for 
the  dispersive  permittivity  and  to  inverse  transform  the  weighted  result. 
The  development  of  a  space-time  integral  equation  can  be  reduced  to  the 
development  of  a  Greens  function  in  the  time  domain  for  a  dispersive  medium. 
To  do  this  a  physical  model  for  the  dispersion  would  be  helpful  and  perhaps 
necessary.  To  date  a  space-time  integral  equation  for  general  dispersive 
media  has  not  been  developed.  The  problem  solution  could  be  carried  out 
by  solving  Maxwell's  equations  directly  or  by  solving  the  frequency  domain 
integral  equation  and  inverse  transforming  the  result. 
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SECTION  III 


SPACE-TIME  INTEGRAL  EQUATION  SOLUTION  -  COMPLEX  CONDUCTING  TARGET 

The  space-time  integral  equation  formulation  is  valid  for  any  target 
shape.  Solution  techniques  have  been  developed  for  a  number  of  important 
shapes  for  the  case  of  conducting  targets  [5j[6].  In  particular,  flat  plates' 
and  cylinders  have  been  treated.  These  shapes  have  also  been  combined  to  form 
the  solution  to  scattering  from  an  aircraft  or  missile  model..  In  this  section 
the  solution  for  a  Tomahawk  Cruise  missile  model  is  given.  The  model  is  a  long 
thin  sphere-capped  cylinder  with  flat  rectangular  plates  for  wings  and  stabil¬ 
izers.  This  shape  very  closely  resembles  the  missile.  The  theory  of  the  solu¬ 
tion  is  reviewed  in  Subsections  3 . 1  and  3.2;  followed  by  the  computed  results 
and  comparison  with  measurement. 

3.1  THEORY  OF  THE  SPACE-TIME  INTEGRAL  EQUATION  SOLUTION  FOR  SCATTERING  FROM 
A  CONDUCTING  MISSILE  MODEL 

The  problem  of  determining  the  scattering  by  cylinders  with  fins  attached 
is  of  great  practical  interest,  since  this  serves  as  a  model  for  numerous  mis¬ 
siles  and  aircraft.  Basically  this  approach  consists  of  developing  two  simul¬ 
taneous  space-time  integrodif ferential  equations  and  their  subsequent  computer 
solution  by  marching  on  in  time.  These  equations  contain  terms  which  may  be 
interpreted  as : 

(a)  The  influence  of  cylinder  currents  on  other  cylinder  currents. 

(b)  The  influence  of  fin  currents  on  cylinder  currents. 

(c)  The  influence  of  fin  currents  on  other  fin  currents. 

(d)  The  influence  of  cylinder  currents  on  fin  currents. 

The  neighborhood  along  the  line  where  a  fin  is  attached  to  the  cylinder  is 
accounted  for  by  application  of  boundary  conditions  at  the  edge  of  the  fin. 

The  boundary  condition  for  a  conductor  is  that  n  x  E  =  0.  In  addition  it 

is  true  that  H  =  0  on  the  surface.  The  equation  for  the  total  field  in  terms 
n 

of  its  surface  values  then  becomes  (compare  (2-1)  ) 


where 


H(r,t) 


-+i 

H  (r,t) 


L{ n’x  t") 


x  R  dS' 


H '  ~  H  ( r  1  ,  T ) 

r  119 

L  ”  2  +  Rc  9l 

K 

*+  -*•  -+ 

R  =  r  -  r' 

T  =  t  -  R/C 


(3-1) 


We  will  normalize  time  to  light-meters  (c  =  1) .  On  the  surface  of  a  conductor, 
(2-11)  simplifies  to 


J(r,t)  = 


-*i  -► 

2  J  (r,t)  + 


2TT 


£ 


L  J'  x  R  dS' 


(3-2) 


The  above  is  the  result  of  applying  a  limiting  process  to  (3-1)  in  moving  the 
observation  point  r  to  S.  We  proceed  as  in  Section  2.  An  incident  field  is  ap¬ 
plied  of  form 


(3-3) 


and  (3-2)  is  first  solved  numerically  for  all  time  by  a  marching-in-time  pro¬ 
cedure.  Having  found  the  surface  currents,  the  farscattered  field  is  obtained 
by  application  of  the  asymptotic  form  of  (3-1) 


(3-4) 


Equation  (3-2)  for  J  on  S  is  applied  on  the  surface  of  a  conducting  solid.  For 
the  thin  conducting  surface  a  different  formulation  is  used. 

The  integrodif ferential  equation  for  the  fin  currents  is  obtained  by 
starting  with  the  expression  for  the  electric  and  magnetic  potentials.  The 
total  electric  field  is  given  by 
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-  V<f>(r,t) 


(3-5) 


F 


I 


i 


4 


-*■  ->  -*-i  -> 

E(r,t)  =  E  (r#t) 


9A(r,t) 

'  y  9t 


where 

A  =  magnetic  vector  potential  such  that  H  =  v  x  A 
$  =  electric  potential 
y  =  permeability  of  space. 

Next,  apply  to  (3-5)  the  Lorentz  gauge  relation 

V  •  a  +  e  =0 


where  £  is  the  permittivity  of  space,  yielding 

£  2EJ|d=L  =  £  +  V(V  •  A)  -  y£  &!f£L  . 

9t 


(3-6) 


Then  applying  the  E-field  boundary  condition  that  the  component  of  the  total 
E-field  tangent  to  the  fin  vanishes  everywhere  on  the  fin,  (3-6)  can  be  ex¬ 
pressed  as 


92A(r,t) 


9t 


9Ei(r,t) 

e 


(3-7) 


or,  in  a  more  simplified  form 


2t 
□  A 


-*• 

A  = 


9E1 

9t 


(3-8) 


where 


i 
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the  wave  operator. 


To  form  the  above  expansion,  the  fins  are  assumed  to  lie  in  the  x  -  z  plane. 

In  order  to  obtain  an  expression  for  the  fin  currents  directly,  it  is  neces- 
sary  to  express  the  magnetic  vector  potential  A  as  the  sum  of  two  components, 
one  due  to  the  observation  point  ("self"  patch)  and  the  other  due  to  the  re¬ 
mainder  of  the  patches.  If  the  fin  grid  patches  are  small  enough,  then  the 
current  over  a  given  patch  can  be  assumed  constant.  This  makes  it  possible  to 
express  A  as 


A(r,t) 


VF£.t) 


/  !“!  -  • 

Non-self 

Patches  T  =  t  -  R 


(3-9) 


where 

-V  ¥  -  (If  -  *)/’• 

As  =  area  of  patch  containing  r 
6  =  thickness  of  fin 


Substituting  this  into  (3-8)  yields  the  fin  currents  directly  as 


(3-10) 


2~*  1  I  2  -+  -*■ 

‘  JF(r'fc)  +  ^'JF(r,t)  =  Y  |  "  E  IT  '  '  ’  V(r,t> 


■+  .  I 

,(r,t)  ) 


where 


.i  / 

4TT  y 


Non-Self 

Patches 


|j(rV 

{  R 


T  =  t  -  R 


The  surface  currents  contributing  to  A,  „  are  both  cylinder  currents  and  fin 

NS 

currents.  The  contribution  of  both  of  these  can  be  separated  in  (3-10)  and 
displayed  explicitly  to  give  the  expression  for  the  current  flowing  on  the  fin 


_ ,2*^  _ _ ,  — ►  — v 

□  Jp(r,t)  +  tag  Jp(r 


'  *7  it '  -2  [\s<J't»  +  V*'*’] 


-  £3  A  (r,t)  + 

FNS 


Ac(?.t,J  J 


(3-11) 


where  Ac(r,t)  is  the  magnetic  vector  potential  due  to  the  cylinder  current. 

The  fin  currents  at  the  free-space  edges  and  at  the  fin-cylinder  join 
are  given  by  the  boundary  conditions 


J  (r,t)  -  0  ,  J„  (rf t)  -*• 


(free-space  edge) 


(3-12) 


9^  (r,t) 


—  Of  >Jf/  —  0 


(free-cylinder  join) 


where  (r,t)  =  current  component  perpendicular  to  the  edge  or  join 

j„  (r,t)  =  current  component  parallel  to  the  edge  or  join. 

This  completes  the  solution  for  a  target  consisting  of  a  solid  with  thin 
plates  attached.  It  should  be  noted  that  while-  separate  equations  are  used  to 
solve  for  on  the  cylinder  and  for  ”3p  on  the  fins,  the  integrals  in  (3-2) 
and  (3-11)  over  ”3  refer  to  both  Jc  and  Jp  so  that  the  equations  are  coupled. 
The  inarching  in  time  procedure  allows  the  solution  by  solving  first  (3-2) 
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and  then  (3-11)  at  each  time  step. 

3.2  NUMERICAL  IMPLEMENTATION 

To  obtain  a  solution  of  the  space-time  integral  equations  for  a  cylinder 
with  fins  attached,  each  of  the  components  of  the  problem  must  be  represented 
mathematically.  First,  the  geometry  of  the  scatterer  and  the  characteristics 
of  the  incident  field  must  be  numerically  specified.  From  these,  the  surface 
and  fin  currents  can  be  computed  using  numerical  representations  of  (3-2)  and 
(3-11) .  This  is  accomplished  by  carrying  out  the  integration  and  differentia¬ 
tion  numerically  and  using  a  "marching  on  in  time"  procedure.  The  far-scattered 
field  can  be  computed  directly  from  the  current  densities  by  using  a  numerical 
representation  of  (3-4) , 

To  describe  the  scatterer  geometrically,  the  scattering  surface  is  div¬ 
ided  into  curvilinear  patches  of  approximately  equal  area  with  a  space  sample 
point  at  the  center  of  each  patch.  The  spacing  of  these  sample  points  (and 
thus,  the  size  of  the  patches)  is  chosen  Small  enough  to  give  both  a  good 
representation  of  the  scatterer  itself  and  of  the  currents  that  exist  on  the 
scatterer.  The  sample  point  spacing  affects  the  time  increment  At  at  which  the 
current  densities  can  be  calculated.  The  time  increment  must  not  exceed  the 
time  it  takes  a  wave,  moving  at  the  speed  of  light,  to  travel  between  the 
closest  space  points.  This  insures  that  the  integral  equation  can  be  expressed 
as  a  recurrence  relation  in  time  and  that  a  matrix  inversion  is  not  necessary 
to  obtain  a  solution. 

The  space-time  integral  equations  represent,  in  principle,  the  solution 
of  the  scattering  problem  for  cylinders  with  fins  attached  for  an  arbitrary 
incident  field.  Although  the  equations  can  be  solved  for  each  incident  field 
separately,  it  is  very  inefficient  to  do  so.  Since  in  most  practical  scatter¬ 
ing  problems  the  excitation  is  a  plane  wave,  a  more  efficient  way  to  approach 
this  problem  is  to  compute  the  scattered  response  when  the  incident  wave  is  an 
electromagnetic  impulse.  Once  the  impulse  response  of  a  target  has  been  ob¬ 
tained,  the  response  due  to  any  incident  plane  wave,  whose  spectrum  is  con¬ 
tained  within  the  spectrum  of  an  impulse,  can  be  calculated  by  a  simple 


convolution  procedure.  Moreover,  the  impulse  response  is  intimately  related  to 
the  actual  geometry  of  the  target,  and  thus,  the  potential  for  developing 
techniques  to  determine  the  impulse  response  of  a  scatterer  by  an  inspection 
of  its  geometry  is  ever  present. 

For  the  numerical  solution,  however,  it  is  not  practical  to  use  an  ideal 
impulse  for  an  excitation.  Thus,  in  this  solution  a  regularized  or  smoothed 
impulse  is  used.  The  form  of  this  illumination  at  the  origin  is  the  Gaussian 
regularization  of  an  impulse,  namely 

H1(t)  =  —  e 
/n 

which  converges  to  a  delta  functional  as  n  goes  to  infinity.  The  time  domain 
integral  equations  can  be  solved  exactly  for  bodies  with  linear  dimensions  up 
to  several  pulse  widths  of  this  regularized  impulse.  In  this  solution  consid¬ 
eration  is  limited  to  bodies  of  this  size. 

The  currents  flowing  on  the  cylinder  and  on  the  fin  are  computed  using 
numerical  representations  of  (3-2)  and  (3-11).  For  the  purpose  of  discussing 
the  numerical  solution,  only  the  x-components  of  these  equations  will  be  con¬ 
sidered.  similar  representations  are  used  for  the  calculation  of  the  other 
components . 


For  the  numerical  solution,  (3-2)  for  the  x-component  of  the  cylinder 
current  at  sample  point  i  and  time  t  is  represented  as 


where 


(1  +  e)  =  self-term  correction  factor  for  observation  patch  i 


(See  Appendix  7-3) 

=  the  y-  and  z-components  of  the  incident  field  at  patch  i, 
time  t 
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r 


N 

F(J) 


J 

x 


T 


t 

R 


n  . 
ni 


Rx 


Ry 


Rz 


number  of  grid  patches  on  the  cylinder  and  fin 
j(r£,x)/R  +  3j (r£>Tj  /3t 

the  x-component  of  the  current  density  at  patch  £,  time  t 
t  -  R 


time  in  light-meters 

distance  from  the  integration  patch  £  to  the  observation 

patch  i:  R  =  /(*.  -  +  (y.  -  y^2  +  (T  -  zzf 

unit  normal  at  patch  i,  In  n  . ,  n  .  \ 

\  xi  yi  zi/ 

unit  vector  from  patch  £  to  patch  i,  /n„  ,  n„  ,  n„ 
v  \  Rx  Ry  Rz 


R 


AS^  =  area  of  patch  £  . 


The  time  differentiation  and  interpolation  necessary  for  the  evaluation 
of  the  integrands  appearing  in  (3-2)  are  performed  numerically  by  represent¬ 
ing  the  surface  current  with  a  fourth-order  polynomial.  In  order  to  achieve 
the  best  accuracy,  the  five  points  used  for  the  representation  are  chosen 
such  that  the  current  is  evaluated  as  near  as  possible  to  the  middle  of  them. 

The  numerical  representation  of  (3-11)  for  the  x-component  of  the  fin 
current  at  ri»t  is  obtained  by  representing  the  second  time  derivative  of 
the  fin  current  by  a  three-point  difference  approximation,  yielding 


(3-15) 


l 
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where 


At  =  solution  time  step 

A  =  y  CT  (r.  ,t  -  At\  +  A 
x  F.  F  \  /  x, 


r  J„(v) 


NS 


NS 


£ 

Z=  1 

VI 


As, 


T  =  t  -  At  -  R  . 


The  space  and  time  derivatives  in  (3-15)  are  evaluated  using  three-  and  five- 
point  difference  approximations.  For  the  time  interpolation  a  linear  approxi¬ 
mation  is  used. 

The  equations  for  the  surface  and  fin  currents  are  solved  with  a  digital 
computer  by  simply  marching  on  in  time.  The  computation  starts  at  a  point  in 
time  before  the  incident  field  reaches  the  scatterer  and  proceeds  sequentially 
in  time  in  the  same  manner  that  nature  would  solve  the  problem  in  the  real 
world.  It  is  important  to  note  that  since  the  minimum  spacing  between  space 
sample  points  on  the  surface  is  not  less  than  At,  then  (3-14)  and  (3-15)  give 
the  current  density  in  terms  of  other  currents  at  times  not  later  than  (t  -  At) 
which  are  already  known.  Thus,  the  integral  equations  (3-2)  and  (3-11)  have 
been  reduced  to  recurrence  formulas  in  time  and  the  need  to  perform  matrix  in¬ 
versions  has  been  eliminated. 

Once  the  current  densities  have  been  computed,  the  far-scattered  field 
can  be  calculated  directly  using  a  numerical  representation  of  (3-13) .  For 
computation  purposes,  (3-13)  is  expanded  in  rectangular  coordinates  yielding 


where 


HS  =  HS  x  +  HS  y  +  HS  z 
x  y  z 

r  =  n  x  +  n  y  +  n  z  . 
rx  ry  rz 

The  time  differentiation  and  interpolation  necessary  for  the  evaluation  of 
3J  r^,T  /3t  are  performed  numerically  using  a  five-point  Lagrange  interpola¬ 
tion  formula. 


The  geometry  parameters  used  for  the  description  of  the  far-scattered 
field  are  illustrated  in  Figure  2-1.  In  this  sketch  the  incident  field  is 
shown  propagating  in  the  yz-plane  and  making  an  angle  a  with  the  z-axis.  The 
scattered  fields  produced  by  the  surface  currents  are  then  computed  in  the  two 
principal  planes.  The  scattered  field  is  computed  in  the  yz-plane  at  angles 

with  respect  to  the  direction  of  propagation  of  the  incident  wave.  In  this 
plane  the  two  orthogonal  components  used  to  represent  the  scattered  field  are 

-►g 

the  component  perpendicular  to  the  yz-plane,  H  ,  and  the  component  tangent  to 
-*s  yzx 

the  yz-plane,  H  The  other  plane  in  which  the  scattered  field  is  computed  is 

the  px-plane,  which  is  formed  by  the  direction  of  propagation  of  the  incident 
wave  and  the  x-axis.  In  this  plane  the  two  components  used  to  represent  the 

g 

scattered  field  are  the  component  perpendicular  to  the  px-plane,  H  ,  and  the 

s  pxp 

component  tangent  to  the  px-plane,  H  These  scattered  fields  are  computed 
at  angles  ,  which  are  measured  with  respect  to  the  direction  of  propagation 
of  the  incident  wave. 


3.3  RESULTS  FOR  MISSILE  MODEL 

A  sketch  of  the  model  used  for  calculation  and  measurement  is  shown  in 
Figure  3-1.  The  dimensions  are  given  in  inches.  The  size  of  the  model  was  cho¬ 
sen  for  convenience  as  L  =  26.25",  D  =  2.5",  this  represents  about  1/8  scale. 
The  relative  dimensions  of  L,  D,  the  wings  and  stabilizers  are  close  to  the 
actual  dimensions  (as  published  in  Jane's  "All  the  World's  Aircraft”).  In  these 
as  in  all  other  respects,  the  model  closely  resembles  the  actual  missile,  ex¬ 
cept  for  some  detail  not  deemed  of  great  importance  to  the  response.  The  actual 
missile  dimensions  are  given  as  L  =  5.5  m  and  D  =  .53  m. 
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The  computed  responses  are  given  in  Figure  3-1',  at  several  aspects  t rom 
a  =  0°  (nose  on)  to  a  =  90c  (broadside).  The  angle  •»  is  measured  in  the  sym¬ 
metry  plane,  which  is  the  vertical  plane  for  a  horizontally  t lying  missile. 

The  time  scale  is  in  units  of  light- inches ,  corresponding  to  the  dimen¬ 
sions  given  for  the  model.  The  features  of  the  response  can  thus  be  related 
directly  to  the  physical  model. 

Consider  the  response  at  a  =  0.  First  a  small  doublet  resembling  the  de¬ 
rivative  of  the  imput  pulse  appears  at  t  =  -25.  This  is  the  specular  response 
from  the  rounded  nose.  The  large  positive  peak  at  t  =  -2.5  corresponds  to  re¬ 
flection  from  the  leading  edge  of  the  wings.  The  remainder  of  the  resjwnse 
results  from  the  complex  interaction  of  wings,  cylinder  and  the  remainder  of 
the  target. 

The  main  effect  in  going  from  a  =  0  to  u  =  90  is  that  the  above  features 
move  closer  in  time.  Finally,  at  a  =  90,  the  specular  response  from  the  cylin¬ 
der  and  that  from  the  wing  surface  join  to  make  a  large  initial  response.  The 
later  bumps  at  t  =  6  and  t  =  12  can  be  interpreted  as  the  responses  due  to 
creep  around  the  fuselage  and  the  wings. 

Measurements  were  made  on  the  time  domain  scattering  range  for  an  alumi¬ 
num  model  of  the  same  dimensions.  The  results  of  these  measurements  are  com¬ 
pared  with  calculations  in  Figure  3-3,  at  several  angles.  It  is  seen  that  the 
comparisons  are  quite  good. 
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FIG.  3-2  Smoothed 


FIG.  3-3  Comparison  of  measured  and  calculated  responses  —  Tomahawk. 


SECTION  IV 


TIME  DOMAIN  MEASUREMENTS 

A  necessary  ingredient  for  any  successful  development  of  techniques  for 
the  solution  of  previously  unsolved  scattering  problems  is  the  availability  of 
actual  measurements  that  can  be  used  to  verify  the  new  techniques.  During  this 
effort  the  Sperry  Research  Center's  Time  Domain  Scattering  Range  was  used  for 
these  measurements.  Section  4.1  describes,  in  summary  form,  this  measurement 
system,  its  operation,  and  its  important  performance  characteristics.  The 
measurement  results  which  were  taken  for  this  effort  are  described  in  Section 
4.2. 


4 . 1  MEASUREMENT  SYSTEM 


The  Sperry  Research  Center's  Time  Domain  Scattering  Range  is  a  system  for 
obtaining  a  low-noise  response  in  the  nanosecond  region.  It  consists  of  a 
ground  plane  scattering  range,  a  step  function  source,  a  sampling  oscilloscope 
receiver,  and  a  laboratory  instrumentation  computer  for  control  and  proces¬ 
sing.  The  system  signal  source  is  a  high-voltage  switch  which  generates  a 
300  V  step  function  with  a  risetime  less  than  100  ps.  The  signal  is  radiated, 
virtually  undistorted,  from  a  vertical  wire  transmitting  antenna  located  at 
the  center  of  a  20  foot  diameter  circular  ground  plane.  This  wave  illuminates 
the  target  and  the  resulting  scattered  waveform  is  received  on  a  flush-mounted 
coaxial  horn  antenna,  which  essentially  smooths  and  differentiates  the  signal 
and  thus  provides  the  smoothed  impulse  response  of  the  target.  The  received 
waveform  is  sampled  by  a  12  GHz  sampling  oscilloscope  that  has  been  triggered 
by  the  initial  pulse  and  whose  sampling  gate  deflection  is  under  the  control 
of  a  small  instrumentation  computer.  Unprocessed  data  are  displayed  on  the 
oscilloscope  CRT  while  the  sampled-and-held  waveform  is  passed  through  a  low- 
pass  filter,  digitized,  read  into  the  computer,  and  stored  on  magnetic  tape 
automatically.  The  waveforms  are  stored  in  such  a  way  that  they  are  ready  for 
the  subsequent  operations  of  averaging  (to  remove  short-term  noise)  and  base¬ 
line  processing.  The  effects  of  baseline  contamination  are  subtracted  from 
measured  waveforms  to  improve  system  accuracy. 
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The  time-domain  scattering  range  provides  a  simple  technique  for  obtain¬ 
ing  transient  time-domain  data  or  multi-octave  frequency-domain  data.  The  cru¬ 
cial  feature  that  the  time-domain  scattering  range  yields,  is  a  "free  time 
window"  between  the  arrival  of  the  direct  wave  and  the  arrival  of  unwanted  re¬ 
flections.  Targets  are  usually  located  anywhere  from  two  to  five  feet  from  the 
transmitting  antenna.  The  response  from  the  antenna  tip  and  the  ground  plane 
edge  occur  at  approximately  15  ns.  Thus,  a  "clear  window"  exists  between  4  ns 
and  15  ns  which  can  be  used  to  view  the  target  responses.  The  entire  region 
between  the  direct  transmission  and  the  range  edge  response  forms  this  "free 
time  window"  for  viewing  the  target  response.  The  unwanted  reflections  are 
gated  out  in  time.  Thus,  undistorted  transient  target  responses  can  be  viewed 
without  resorting  to  anechoic  chambers. 

The  accuracy  of  the  measurement  system  has  been  estimated  for  the  results 
presented  in  this  section.  The  peak  of  the  incident  pulse  as  measured  on  the 
sampling  oscilloscope  is  approximately  400  mV,  and  a  typical  target  response 
has  a  peak  value  in  the  vicinity  of  5  mV.  When  using  the  10  mV  scale  on  the 
sampling  oscilloscope,  the  standard  deviation  of  the  sample  mean  is  estimated 
to  be 

0_  =  0.5  mV 
v 

if  16  scans  are  averaged.  Thus,  the  estimated  standard  deviation  of  the  sample 
mean  v  is  in  the  vicinity  of  10%  of  the  peak  value  of  the  target  response. 

In  addition,  the  measured  responses  are  further  processed  by  means  of  a 
convolution  procedure  to  obtain  the  response  due  to  a  Gaussian  shaped  incident 
pulse  rather  than  the  approximate  smoothed  impulse  used  in  the  actual  measure¬ 
ments.  Figure  4-1  displays  the  actual  measured  incident  pulse  along  with  the 
smoothed  Gaussian  pulse  that  was  used  in  the  convolution  process.  The  frequen¬ 
cy  spectrum  of  these  two  pulses  is  displayed  in  Figure  4-2.  The  time-domain 
expression  for  the  Gaussian  pulse  is  given  by 


and  the  frequency  domain  expression  is 
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E(w)  =  Eq  exp  £-  oa2/|2an|2j  exp^-j(jJtQ 

The  amplitude  coefficient  Eo  was  set  equal  to  the  dc  value  of  the  measured  in¬ 
cident  pulse.  The  width  coefficient  a^  was  obtained  by  requiring  that  the  50% 
value  of  |e(w)|  occur  at  the  same  point  in  frequency  as  the  50%  value  of  the 
magnitude  of  the  transformed  measured  incident  pulse.  The  width  of  the  result¬ 
ing  Gaussian  pulse  becomes  0.59  ns  or  approximately  7  inches.  An  added  benefit 
of  this  process  is  the  reduction  of  high  frequency  noise  in  the  response. 

This  measurement  system  provides  the  capability  of  measuring  and  record¬ 
ing  the  processed  smoothed  impulse  response  data  on  computer  computable  tape. 
These  data  may  then  be  used  to  find  the  response  due  to  any  radar  waveform 
whose  spectrum  is  contained  within  the  spectrum  of  the  original  measurement 
smoothed  impulse  illumination.  Futher  details  of  this  system  may  be  found  in 
the  references  [5,6]. 

4.2  MEASUREMENT  RESULTS 


In  previous  studies  the  smoothed  impulse  response  was  measured  for 
numerous  target  geometries.  These  previous  measurements  are  summarized  in 
reference  [6].  During  the  present  contract,  smoothed  impulse  responses  were 
measured  for  a  dielectric  sphere,  a  dielectric  sphere-capped  cylinder  and  a 
dielectric  right  circular  cylinder.  The  dielectric  material  used  for  the  tar¬ 
get  models  was  nylon  which  has  a  relative  permittivity  of  about  3  [8].  The 
smoothed  impulse  responses  of  the  same  shaped  conducting  targets  were  also 
measured  for  comparison  purposes.  Finally,  the  smoothed  impulse  response  of 
the  Tomahawk  missile  model  displayed  in  Figure  3-1  was  measured.  A  summary 
of  the  target  geometries  used  for  the  measurements  is  displayed  in  Table  I. 
The  smoothed  impulse  responses  which  are  shown  in  the  remainder  of  this 
section  are  due  to  the  incident  smoothed  impulse  shown  in  Figure  4-1,  which 
has  a  width  of  °-59  ns  and  spectrum  with  components  out  to  4  GHz  as 

shown  in  Figure  4-2.  In  viewing  the  measured  smoothed  impulse  response  re¬ 
sults  that  are  displayed  in  Figures  4-3  through  4-7  in  the  remainder  of  this 
section,  it  should  be  kept  in  mind  that  two  standard  deviations  represent 


Spectrum  of  incident  pulse. 
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about  1  mV.  Thus,  apparent  features  of  the  waveforms  that  have  variations 
less  than  this  have  less  statistical  significance  and  should  be  discounted. 

In  Figure  4-3  the  measured  smoothed  impulse  response  of  spherical  targets 
is  displayed.  Both  nylon  and  conducting  spheres  of  8"  diameter  and  4"  diameter 
were  measured.  The  first  positive  pulse,  appearing  slightly  before  1  ns,  is 
due  to  the  specular  return.  This  amplitude  for  the  4"  diameter  spheres  is  ap¬ 
proximately  half  that  of  the  8"  diameter  spheres.  Note  also  that  this  ampli¬ 
tude  of  the  dielectric  spheres  is  approximately  one-third  that  of  the  conduct¬ 
ing  spheres.  The  creeping  wave  return  is  evident  at  t  =  2.5  ns  for  the  8"  con¬ 
ducting  sphere  and  at  t  =  1.8  for  the  4"  conducting  sphere.  The  returns  from 
waves  traveling  through  the  dielectric  spheres  and  being  reflected  from  the 
back  side  are  apparent  at  t  -  3  ns  for  the  8"  dielectric  sphere  and  at 
t  =  2.1  ns  for  the  4"  dielectric  sphere.  The  presence  of  baseline  contamina¬ 
tion  of  1  mV  and  less  also  appears  to  be  present  at  later  points  in  time. 

The  measured  smoothed  impulse  response  of  a  nylon  and  a  conducting 
sphere-capped  cylinder  with  TE  polarization  are  shown  in  Figure  4-4.  The  cyl¬ 
inder  has  a  length  of  12  inches  and  a  diameter  of  4  inches.  Axial  incidence 
corresponds  to  zero  degrees.  The  target  is  rotated  about  its  center  for  these 
measurements  which  correspond  to  the  time  t  =  3  ns  on  this  figure.  The  conduc¬ 
tor  responses  in  Figure  4-4 (b)  are  included  for  completeness  where  the  specu¬ 
lar  return,  small  side  return,  and  creeping  wave  return  can  be  easily  discer¬ 
ned.  The  dielectric  measurements  displayed  in  Figure  4-4 (a)  are  of  low  ampli¬ 
tude  and  in  most  cases  less  than  1  mV.  Again,  features  similar  to  those  of  the 
conducting  case  may  be  seen  in  Figure  4-4 (a)  but  other  variations,  especially 
in  the  0°  and  30°  results  later  in  time,  may  be  due  to  some  systems  error. 

Figure  4-5  displays  the  measured  smoothed  impulse  response  for  both  di¬ 
electric  and  conducting  right  circular  cylinders  with  TE  incidence.  The  cyl¬ 
inder  has  a  length  of  8  inches  and  a  diameter  of  12  inches.  Axial  incidence 
corresponds  to  zero  degrees.  The  target  is  rotated  about  its  center  for  these 
measurements  which  correspond  to  the  time  t  =  2.5  ns  on  this  figure.  Both  the 
dielectric  and  conducting  results  are  of  good  quality.  Again  the  specular  re¬ 
turn  amplitude  for  the  dielectric  is  approximately  one-third  the  value  meas¬ 
ured  for  the  conducting  target.  The  portion  of  the  dielectric  return  due  to 
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FIG.  4-3  Smoothed  impulse  response  of  spheres  (hor.  scale:  1.0  ns/div.; 
vert,  scale:  5mV/div.). 


4-4  Smoothed  impulse  response  of  sphere  capped  cylinder  with  TE  polarization 
(hor.  scale:  1.0  ns/div.;  vert,  scale:  5  mV/div.). 
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FIG.  4-5  Smoothed  impulse  response  of  right  circular  cylinder  with  TE  polarization 
(hor.  scale:  1.0  ns/div.;  vert,  scale.  5  mV/div.). 


(a)  NYLON  SPHERE  CAPPED  CYLINDER  (c)  NYLON  RIGHT  CIRCULAR  CYLINDER 


(b)  CONDUCTOR  SPHERE  CAPPED  CYLINDER 


(d)  CONDUCTOR  RIGHT  CIRCULAR  CYLINDER 
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FIG.  4-6  Smoothed  impulse  response  of  cylinders  with  TM  polarization 
(hor.  scale:  1.0  ns/div.;  vert,  scale:  5  mV/div.). 
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wave  traveling  through  and  being  reflected  is  more  complicated  than  the  creep¬ 
ing  wave  effects  noted  in  the  conductor.  In  fact,  these  returns  in  the  dielec¬ 
tric  case  are  probably  a  number  of  waves  which  have  been  called  glory  waves. 

The  TM  responses  for  broadside  incidence  for  both  the  sphere-capped  cyl¬ 
inder  and  the  right  circular  cylinder  are  shown  in  Figure  4-6.  As  before,  the 
specular  return  for  the  dielectric  target  is  approximately  one-third  of  that 
for  the  conducting  target.  In  this  case,  however,  the  internally  reflected 
wave  of  the  dielectric  is  about  twice  the  size  of  the  creeping  wave  of  the 
conductor.  It  also  has  a  more  complicated  structure. 

Finally,  Figure  4-7  displays  the  smoothed  impulse  response  of  the 
Tomahawk  missile  model  that  is  sketched  in  Figure  3-1.  Nose-on  incidence  cor¬ 
responds  to  zero  degrees.  The  target  is  rotated  about  its  center  which  cor¬ 
responds  to  time  t  =  5  ns  on  this  figure.  Note  also  that  the  time  window  is 
10  ns  here  rather  than  the  5  ns  in  the  previous  results.  The  primary  returns 
from  the  nose,  wings  and  tail  region  can  be  easily  distinguished.  The  second¬ 
ary  returns  from  the  traveling  currents  and  interactions  can  also  be  observed. 


FIG.  4-7  Smoothed  impulse  response  of  Tomahawk  missile  model 
(hor.  scale:  1.0  ns/div.;  vert,  scale:  5  mV/div.). 
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TABLE  I 

SUMMARY  OF  TARGET  GEOMETRIES 

8"  diameter  sphere 

(1)  nylon  (er  =  3) 

(2)  conducting 

4"  diameter  sphere 

(1)  nylon  (£r  =  3) 

(2)  conducting 

4"  diameter  by  12"  long  sphere-capped  cylinder 

(1)  nylon  =  3) 

(2)  conducting 

4"  diameter  by  8"  long  right  circular  cylinder 

(1)  nylon  (er  =  3] 

(2)  conducting 

Tomahawk  missile  model 
(see  Figure  3-1) 


SECTION  V 


SUMMARY  AND  CONCLUSIONS 


A  time  domain  solution  was  formulated  and  implemented  for  the  problem 
of  electromagnetic  scattering  from  a  dielectric  solid.  The  approach,  called 
the  space-time  integral  equation  approach,  is  applicable  to  bodies  of  any 
shape.  The  computed  time  domain  response  can  be  directly  related  to  target 
shape.  In  the  frequency  domain  the  equivalent  range  of  validity  is  for  ka 
from  0  to  6  or  more  (well  into  the  resonance  region) ;  techniques  to  extend 
this  over  the  entire  spectrum  have  been  presented  for  conducting  targets  in 
[ 5 1 -  The  approach  was  demonstrated  on  a  dielectric  sphere  and  sphere  capped 
cylinder.  The  solution  for  the  dielectric  sphere  was  verified  by  comparison 
with  the  known  classical  solution.  The  solution  for  the  sphere-capped  cylin¬ 
der  was  verified  by  measurement.  The  problem  of  extending  the  solution  to 
targets  of  dispersive  materials  was  discussed. 

The  technique  was  applied  to  a  conducting  target  of  complex  shape 
resembling  closely  a  Tomahawk  missile.  Responses  were  calculated  and  veri¬ 
fied  by  measurement.  The  time  domain  responses  were  shown  to  be  highly 
suggestive  of  principal  target  features. 

The  tine  domain  scattering  range  was  described.  Results  were  given 
for  several  conducting  and  dielectric  targets  at  several  aspects. 
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SECTION  VII 
APPENDIX 


7.1  TOO  EQUATION  STIE  FORMULATION 


A  more  economical  formulation  of  the  dielectric  problem  solution  is  one 
which  utilizes  only  the  H-field  equations  for  the  inside  and  outside,  namely 
equations  2-11,  2-12  and  2-15,  2-16.  Separating  out  the  self  integrations, 
we  obtain  as  in  Section  2.2: 


(?  -  V) 

(r  +  r) 


-y  a.  9m  -+i  ~  -y 

J  +  YnxaT=2J  +  nxI 


-y  a.  9M  ~ 

J  "  Y  £r  n  X  3T  =  _n  X  Z1 


(7-1) 


For  H  we  can  take  the  average  of  2-15  and  2-16: 
n 


2  H  +  n 
n 


(7-2) 


The  system  (7-1)  can  be  solved  for  J  and  3M/3t.  For  surface  components  Ju 

and  J  : 
v 


where  e  =  y  (K  -  k  )/4. 

u  v 

This  is  just  the  same  equation  as  the  first  of  (2-34)  with  6  =  l/er- 
The  difference  is  that  here  we  expect  to  obtain  3M/3t  by  back  substitution 
in  (7-1) .  There  is  difficulty  in  this  formulation  in  that  the  role  of 
3M/3t  in  the  equations  is  proportional  to  the  patch  radius  y. 


The  results  of  encoding  this  approach  were  poor:  for  some  values  of 
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£  obviously  incorrect  results  were  obtained,  while  for  others  the  solution 
was  unstable.  It  is  not  apparent  why  this  was  so.  Since  the  economy  of  this 
approach  was  not  very  great  it  was  abandoned  in  favor  of  that  of  Section  2. 

7.2  PLANE  SYMMETRY  CONDITIONS 

The  targets  that  have  been  considered  all  have  at  least  plane  symmetry. 
Therefore  the  program  was  coded  with  the  assumption  of  symmetry  about  the 
y-z  plane.  This  reduces  the  amount  of  storage  required  by  a  factor  of  2 
and  also  reduces  the  computation  time  by  a  factor  of  2.  The  incident  direc¬ 
tion  is  always  taken  in  the  plane  of  symmetry  (see  Fig.  2-1)  .  The  following 
table  summarizes  the  symmetry  conditions  on  the  target  surface.  A  plus  sign 
means  that  the  component  sign  is  not  changed  upon  reflection  in  the  y-z 
plane;  a  minus  sign  that  it  does  change;  and  a  zero  that  the  component  is 
zero. 


TE 

Case 

TiM  Case 

Quantity 

X 

y  z 

x  y  z 

-> 

r 

- 

+  + 

-  +  + 

0 

+  + 

+  00 

E1 

+ 

0  0 

0  +  + 

n 

- 

+  + 

+  + 

■y 

J 

+ 

-  - 

+  + 

-y 

M 

— 

+  + 

+ 

H 

scalar  + 

scalar  - 

n 

E 

scalar  - 

scalar  + 

n 

7.3 

SELF  TERM  INTEGRATIONS 

The  numerical  implementation 

of 

the  integral  equation  solution  (equa- 

tions 

;  2-11  through  2-18  of  Section 

2. 

.  1)  is 

made 

more  accurate  by  the  inclu 

sion 

of  terms  obtained  by  analytic 

integration 

on  the  self-patch.  We 

rederive  the  results  of  these  integrations 

here 

.  We  will  evaluate  terms 
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of  the  form 


Since  M  lies  in  the  surface,  the  quantity  n  •  M(r',t)  is  odd  about  r.  Hence 


N1S  =  0  • 


The  remaining  self  integrals  require  some  relations  from  differential  geometry. 


First  consider  the  scalar  integral 


=  S  .  -L  /  /-L  +  h 

2’  4  Ir2  R3t;  ” 


R  dS'  . 


We  need  an  expression  for  (n  ■  R)  ,  where  R  =  (r -r'j/R.  Choose  a  principal 
coordinate  system  u,  v,-  so  that  the  Cartesian  components  of  a  point  on  the 
surface  are  r  =  (x(u,v),  y(u,v),  z(u,v)).  At  observation  point  r  define 


ru  3u  '  rv  3v  '  ‘uu  .  2  '  ''w  .  2 


y2-+  .2-*- 
3  r  -+  3  r 


(For  a  principle  coordinate  system  we  have  that  r  •  r  =  0  and  r  •  n  =  0.) 
„  J  u  v  uv 

r  and  r  are  surface  vectors,  with  r  •  n  =  0  and  r  •  n  =  0.  The  principal 

u  v  u  v  e 


curvatures  are  defined  by 


K  =  (r  ’  n\/(r  •  r  \ 

u  |  uu  Mu  u  / 

K  =  (r  •  n\//r  •  r  j 

v  \vv  f  \  v  v I 


R  can  be  expanded  as 


-v  ->  2.  2  1  ->  2  -* 

-R  =  r 1  -  r  =  r  du  +  r  dv  +  —  r  du  +  —  r  dv  +  r  du  dv 
u  v  2  uu  2  w  uv 


Thus  we  have  that 


(n  •  R)  = 


2R  (KJ*J 


2  2  |  1 2  2 
du  +  K  r  I  dv 
v  v 


W»i  1 1  I 


Expand 


H '  =  H  (r\*r)  -  H  (r ,  t)  +  R  — ~  +  (t  -  t) 

n  n  n  9  R  dt 


and  assume  that  9H  /9R  =  0  and  9H  /9t  (r,t)  =  9H  /9t  (r',T);  then 
n  n  n 


(r2  +  R3t)Hn  (?,t) 


Taking  the  constant  quantities  outside  the  integral: 


(k  I?  I2  du2  +  K  Ir  I2  dv2) 


2S  4tt  n 


Integrating  over  a  circular  patch  of  radius  y  =  /  AS/tt  ,  letting 

Ir  I  du  =  R  cos  0,  Ir  |  dv  =  R  sin  0,  dS '  =  R  dR  d0  t  we  obtain 
1  u  v 

N  =  -  ~~  H  -f-  |  K  cos2  0  +  K  sin2  0j  dR  d0 
2S  4tt  n  /  I  u  v  / 

A  c  '  ' 
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Since  f  cos  0  d0  =  i, 
0 


/K  +  K  \ 

*  -n  4^)  • 


We  have  used  elsewhere  the  definition 


/K  +  K  \ 

-A-^) 


The  vector  version  of  this  integral  is 


R  dS’ 


Here  we  have  to  form  n  x  R.  We  see  readily  that  n  x  R  is  an  odd  function 
— ► 

of  r ' .  Hence 


>r fn.r- 


The  last  of  these  self  integrals  is 


-V 


I 


3S 


x  R  dS'  . 


Write  J'  in  u,v  components 


J'  =  j*  a  +  J'  a  , 
u  u  v  v 


where  a  =r/|r|,a  =r/|r|.  Then,  using  (7-9): 

u  u  u1  v  v  1  v1 


n  x  /  J'  a  x  r\  =  J '  — ^  (k  I  r  |2du 
\  u  u  I  u  2R  y  u1  u1 


2  -  K  |r  |2  dv2 

V1  V1 


and 


n  x  ( J '  a  xr)=J'-|(k  |r  |2  dv2  -  K  | r  | 2  du2) 

\  V  V  I  v  2R  ^  v'  V1  u1  u1  / 

Again,  we  make  the  same  assumptions  as  in  (7-10)  for  the  constancy  of  ^ 
namely  that  3Ju/3R  =  0  and  3Ju/3t  =  constant,  so  that 


( j_  +  jlU.  „  Mi 

W  R3V  r2 


t) 


The  integration  is  similar,  except  that  now  we  get  a  factor  (K  -  IO , 
moreover,  this  factor  is  of  opposite  sign  for  the  a^  and  a^  components 
The  result  is 


<r,t) 


e ("a  J  (r,t) 
u  u 


a 

v 


Jv (r ,t)J 


(7-15) 


(7-16) 

(7-17) 

t 

(7-18) 

and 

(7-19) 


where 


A  "T  A 

The  scalar  version  of  this  integral  is  zero  because  n  •  J  x  R  is  an  odd 
function: 


N  =  n 


1 


f  (-hr+  rlrA  J’  x  R  dS1  =  0  . 


MISSION 

of 

Rome  Air  Development  Center 

RAVC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  In  support  o£  Command,  Control 
Communications  and  Intelligence  (C3I)  activities.  Technical 
and  engineering  support  within  areas  ob  technical  competence 
Is  provided  to  ESV  Program  Ob  bices  IPOs)  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  ob  ground  and  aerospace  objects.  Intelligence  data 
collection  and  handling,  Inbormatlon  system  technology. 
Ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


